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Abstract. We investigate the properties of an axisymmctric non-magnctizcd gas flow without 
angular momentum on a small compact object, in particular, on a Schwarzschild black hole in 
the supersonic region near the object; the velocity of the object itself is assumed to be low 
• compared to the speed of sound at infinity. First of all, we see that the streamlines intersect 

I (i.e., a caustic forms) on the symmetry axis at a certain distance from the center on the 

■ front side if the pressure gradient is neglected. The characteristic radial size of the region, 

in which the streamlines emerging from the sonic surface at an angle no larger than Oq to 
the axis intersect, is Ar = r^O^/S. To refine the flow structure in this region, we numerically 
compute the system in the adiabatic approximation without ignoring the pressure. We estimate 
the parameters of the inferred region with anomalously high matter temperature and density 
' accompanied by anomalously high energy release. 

CN ■ 

<N ; 1. INTRODUCTION 

^ ■ Bondi-Hoyle accretion is the fall of matter onto a moving compact object. Although this 

. problem was formulated in the mid-20th century (Bondi and Hoyle 1944), many details of 
this process are still incomprehensible. In particular, there is virtually no detailed information 
I about the properties of a flow in the supersonic region near a gravitating center in the model 
O ■ of smooth passage of the sonic surface (Bondi 1952). The first approximation to spherically 
^ . symmetric accretion was first considered to calculate the accretion onto a nonrotating compact 
object moving through a gas in the article (Beskin and Pidoprygora 1995). In this case, the 
P [ ratio of the object's velocity foo to the speed of sound at an infinite distance from the body Cqo 
plays the role of a small parameter: 



t3 



We consider only the first approximation with e < 1. This inequality holds for certain astro- 
physical objects. Accretion without shock formation is possible in this case. Since spherically 
symmetric smooth transonic accretion was shown (Bondi 1952; Garlick 1979) to be stable, it 
has a physical meaning. It can then be assumed that Bondi-Hoyle accretion without shock 
formation is also stable at fairly small e and, hence, also has physical meaning. 

Recall the main properties of the smooth spherically symmetric solution denoted by the 
superscript (0) and the first approximation to it denoted by the superscript (1) (Bondi 1952; 
Beskin and Pidoprygora 1995). The solution is sought in the adiabatic approximation with a 
constant adiabatic index P. All of the coefficients k with different subscripts that appear below 
were calculated by Beskin and Pidoprygora (1995) and depend only on the adiabatic index P. 

(1) The requirement of smoothness, i.e., the absence of shocks, leads to an additional condi- 
tion on the sonic surface that yields r^, = ri^^ + ri^^ for the radius of the sonic surface in the 
nonrelativistic limit, where in polar (r, 6, (p) coordinates 

= ^, r« = . . rf fc,(P) cos^. 



1 



(2) The Grad-Shafranov equation defining the stream function $(r, 9) has a solution $ = 
$(0) _^ $(1)^ where 

$(0) = - cos^), = ^o£g(r)sm^e. (1) 

Here 2$o is the accretion rate, the stream function is defined as nv = V$ x e0/27rrsin^, and 
e is a unit vector along the axis. 

(3) The radial function g{r) behaves asymptotically as g{r) = K(T) (r/r*)^ far from the sonic 
surface, which corresponds to a homogeneous flow. On the other hand, in the supersonic region 
r < r* 

y(r) = Mr) (r/r.)-^/% (2) 

where is a radius of the sonic surface. 

(4) We see that the perturbation becomes significant in the supersonic region, and 
max5)($(^Y^o) reaches unity at 

Tx^nieh-af. (3) 

So, the perturbation theory does not hold there and we have to solve exact equations. 

In this paper, we describe the method used and the simplifications that help to carry out 
calculations. We prove that these simplifications can be introduced without any significant loss 
of the accuracy and give basic formulas and results. Subsequently, we discuss whether this 
phenomenon is observable. 

2. DESCRIPTION OF THE METHOD 

The essence of the described method is to directly calculate the streamline, which allows the 
physically observable quantities to be easily found. This method is also the most natural for 
deriving clear equations. 

Let us parameterize the flow line as 9 — 9{9o,r). The initial angle 9o at some radius ro and 
radius r are independent variables. We choose Tq in the supersonic region where it makes sense 
to consider the flrst approximation. The thermodynamic potentials are functions of the same 
arguments. 

First of all, note that this parametrization is unique. However, to specify the flow, we must 
specify not only the trajectory, but also the velocity along it v[9{9o,r),r]. In this case, two 
differential equations for the two functions v and 9 can be derived. They should then be solved. 
These are the energy equation and the equation of forces or angular momcntums. 

The problem is solved by assuming the absence of energy release at the stellar boundary; 
i.e., the object is actually assumed to be a black hole. There is no angular momentum along 
the ^ = axis. 

3. SIMPLIFYING ASSUMPTIONS 

Let us introduce several simplifications none of which, as we will see, affects significantly the 
result. 

(1) In our calculations, we use the metric of flat space. 

(2) The radial velocity of the matter Vr is much higher than the nonradial velocity vq. There- 
fore, the tangential components of the gradients for all parameters of the system are much 
smaller than their radial components. 

(3) To determine the radial velocity, we neglect the enthalpy of the matter compared to the 
gravitational energy. 



(4) The solution can be represented as a converging series in 9q in the form 

oo 
n=0 

on the front side of a compact object near the symmetry axis (6'(0, r) = on the axis). 

4. BASIC EQUATIONS 

Let us first solve the problem by disregarding the pressure in the supersonic region. We go 
to the frame of reference, in which the body is at rest and the gas flows on it. Let us write the 
angular momentum conservation equation for the gas relatively to the center of the object as 

— rV = 0. (4 
dt\dtj ^ ^ 

Next, we can write the expression for the radial velocity as 

dr fr. 



Vfun ~ \vr\ = ~ y 

where Tg = 2GM/c^ is the gravitational radius of an object of mass M. Ehminating dt from 
equations (4,5), we immediately determine the streamline as 

d_ / 3/2^\ _ Q 

dr \ dr / 



The solution of this equation that satisfies the natural initial condition ^(^o, '"o) = is 



e = eo- c{eo) 



roy/^ _ 
r J 



(6) 



We derive the form of the function C{9o) from (2), which yields the same solution for the 
streamline, but has a limited validity range as the first approximation. As a result, we obtain 



6{eo,r) = 6*0 - eki^y/n ( ]=] sin 6*0. 



This expression has a nontrivial structure. It defines a caustic near the distance Rbaii from 
the center, 

-Rball — '"*(£Mn)^- 

For physically meaningful F (F > 1.25), this caustic is located on the front side of the object, 
although the saddle point with a zero gas velocity is always on the rear side. The characteristic 
radial size of the region, in which the streamlines emerging from the sonic surface at an angle 
no larger than Oq to the axis, intersect is Ar = rxO^/S. Naturally, the gas pressure shouldn't be 
ignored in a region of significant compression; the trajectories will not intersect if the pressure is 
included. We can qualitatively predict that, in this case, a certain region of strong compression 
will be closer to the object than i^baii- Next, let us perform a calculation without neglecting the 
pressure. Let us set up a continuity equation for the gas flow. The area of the cross-sectional 
element of the flow in a direction perpendicular to the radial vector is 

5S = 27rr'^ sin 956. (7) 

The continuity equation itself from which we find the number density n can then be written as 



VrTiSS = const. 



(8) 



Using (5) and (7), we obtain from (8) 



sin 656 = ro^'^^no sin 6q56q, 



where 6 is the function of 6'o at Tq. We then get 

Hence, the number density can be written as 

ToX^/^sin^o 1 



n = riQ 



/roy/^sinyo_L (g^ 
\rJ sine doJ' ^ ' 



Below, for simphcity, we consider an adiabatic process on an ideal gas, p — k{s)n^ , where the 
entropy s is constant in the entire space. The pressure, the temperature, and the speed of 
sound can then be expressed using (8). 

We can now calculate the curvature of the streamline under pressure. The point probe mass 
is affected by the gravity force and the force proportional to Vp. We can derive an equation for 
gas angular momentum deviation relatively the body center; the gravity torque is zero indeed. 
The angular momentum is the vector perpendicular to the plane passing through the symmetry 
axis. Only the tangential part of Vp will appear in what follows. 

The moment equation is 

= -AK, 

dt 

where L is the angular momentum, and K is the moment of force. The minus corresponds to 
the repulsion of streamlines. 

The force acting upon the small element is AF = lArAp, where / is the length of the 
circumference described in the plane perpendicular to the ^ = axis. The corresponding 
moment of force is AK = IrArAp, and the mass of the element is Am = mpfilrAeAr, where 
rrip is the mass of a single particle. The angular momentum of this probe mass is AL — 
m-pnlrAOAr ■ r^dO/dt. Finally, we find that 

d f ^2dO\ 



dt \ dt J AOm^n 



At r = const we obtain p = p{9o) and 9 = 6{6q); therefore Ap/ AO = dgQp{9o,r)/dQg9{9o,r). We 
replace the differentiation with respect to t by the differentiation with respect to r using (5): 

\ / mpndg^9{9o,r)rg 



Substitution in the number density in form (9) yields 

2-r 



dr[r'/''dr9{9o,r)] 



sin 0(00 , r) 
sin 00 



-do. 



sin^n 



(11) 

This expression can be rewritten in convenient dimensionless coordinates. Let us denote 



4r^ 



So = 1, and Xq = (ro/r*). As a result, equation (11) takes the form 



|(4_r) 



-da 



sin 6q 



xoj {de^e{eQ,x)) °\sme{eo,x) dg^e{eo,x) 

(11') 

We have obtained a second order partial differential equation of the hyperbolic type that 
is linear in second partial derivatives. Let us transform it to a system of ordinary differential 
equations. For this purpose, we use assumption (4) and substitute 9 with a finite sum containing 
the terms with the powers of 9 up to 2A'" + 1 
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^(^o,r) = J]^o'"+^W(r). (12) 



n=0 



We know the function 9{9Q,r) (10) at the radius Tq which defines initial conditions for (18). 
Then 

A;i(ro) = l, WM = 0, n = 1, 2, A^. (13) 



- (-^'"ap^rfTT^V^- (14) 

The Cauchy problem is Eq. (11) with the initial conditions (13) and (14). After substituting 
(12), we obtain a system for + 1 functions of the radius. The absence of terms with even 
degrees can be explained by the zero initial conditions for them, while the equation to be solved 
is homogeneous. 

The seeming complexity of the method of solution described above can be easily explained 
when it is considered that we must separate out a region of a very small size in both ^ ^ 1 
and r when solving the second-order partial differential equation. Qualitatively, the behavior 
of the trajectories can be described as follows. A plane uniform flow at inflnity transforms into 
an almost spherically symmetric flow near the sonic surface. Subsequently, depending on the 
adiabatic index, three cases are possible (Beskin and Pidoprygora 1995): 

1) . For r — > 5/3 the Mach number M, which deflnes the ratio of the gravitational forces 
to the pressure force, does not tend to inflnity, but retains a value of the order of unity as 
r approaches zero. It can be no signiflcant additional compression compared to that in the 
spherically symmetric case. 

2) . For r < 1.25 the parameter A;in(r) is negative with positive K{T), therefore, the stream- 
lines converge on the rear side (this can be seen from (1)). 

3) . For r > 1.25 the sign of the coefficient k^^ is the opposite; the streamlines converge on the 
front side of the object. After the passage of the region of maximum streamline convergence 
near the symmetry axis, the pressure pushes them apart. 

Thus, our region is described by the following quantities: the radius r^,, at which there 
is maximum additional compression on the axis compared to that in Bondi accretion; the 
dimensionless radius Xx', the radial size of the region Ar at the boundary of which the additional 
compression decreases by a factor of 2; the caustic radius at the same parameters, but without 
including the pressure -Rbaih the minimum achievable ratio 9/9q in this region k^] the minimum 
Mach number in the vicinity of r^- 
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Let us make a calculation for two cases: F = 4/3, — 0.0002 and F = 7/5, = 0.00001. 
The first case is nonphysical. It has only a theoretical significance as an extension of classical 
astrophysical problems. In turn, the sonic surface with the parameters corresponding to the 
second case can be passed in dense clouds of molecular hydrogen. The characteristic radial 
size of the region is Ar ^ rx/5. The solution is stable for r > Tj,; i.e., the result does not 
change with increasing N. Further out, the algorithm is unstable; therefore, the behavior of 
the solution cannot be determined at small r ^ r^,. 

Let k{r) = ki{r), accordingly, = k{r-x)- Let us estimate the physical parameters of the 
system on the axis. First, the number density can be written as 

/To \ 

n(0, r) = no y—j k{r)~ , 

which yields 

1 n^^-'^)/^ 



M = 



^(3r-5)/4^r-i 



and 



for the Mach number. Thus M(0, r) = (rkr/r)^^~^^^^'^k{rf~^ 

(5-3r)/4 

(Fk- l(3r-5)/2, r-1 



-Rball \ 

rx J 



For F = 4/3 we have = 0.025, and for F = 7/5 we have = 0.026 (Beskin and 
Pidoprygora 1995). In fact, the small parameter for ^^^^ is s/cin rather than e, which has an 
order of 0.1 ^ 1. 

Figure 1 shows the dependence of the ratio of the temperatures in our case and in the 
spherically symmetric case on the dimensionless distance to a compact object calculated for 
the above parameters. 

Figure 2 shows the corresponding bivariate dependence of the temperature on the dimension- 
less distance and the angle 9o. The transformation to the {r/vg, 9) plane near is to be made 
by a factor of 1/kx compression towards the symmetry axis of the system. The light spots in 
the lower corners of the figure have no physical meaning, but are determined by the divergence 
of the algorithm. Using Fig. 1, we can associate the brightness in Fig. 2 with temperature. 



5. VALIDITY OF THE ASSUMPTIONS 

To solve the problem, we made many assumptions. Let us substantiate them. For our 
estimates, we set Coo = 0.0002. 

(1) In most cases, the needed region is located far from the event horizon of the object, 
r ^ Tg, which allows us to consider the metric fiat. 

(2) The inclination of the streamline to the radial vector is near r^. The solution for n(r) 
yields Ar = 5 - lO^cm, = 10'' cm, and k^ = 0.03. The region is compressed almost uniformly 
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Figure 1. Temperature ratio along the 9 = axis versus r/vg near the singu- 
larity for the parameters F = 7/5, e = 0.03, = 0.00001; Tq corresponds to 
Bondi (1952) accretion. 

for 00 < 0.3. whence, we can estimate a as a ^ OokxTx/Ar ~ 0.2, and this is the maximum 
possible value for a within our range of parameters. 

(3) To substantiate the possibility of ignoring the enthalpy of matter compared to the grav- 
itational energy, we calculated the Mach number M. The ratio of the enthalpy to the kinetic 
energy for a non-relativistic gas is 2c^/f ^(F — 1) = 2/M^(F — 1), reaching unity at r^. However, 
this is true in a small region. This effect can be easily estimated by solving the equations for 
the effectively increased coefficient of the term responsible for the pressure. When the coeffi- 
cients increases by a factor of 1.5, decreases by approximately the same factor 1.5, while 
increases by no more than 10%. 

4. The smallness of ^o^^^25('"x) compared to Ooki(rx) corresponds to 6o < 0.3 at r = r^:, which 
corresponds to an area fraction of 0.036 of the initial one at tq. This implies that much 
(> 1/30) of the flow is subjected to strong compression. At a large radius, i.e., larger than r^,, 
each succeeding term ^^""'"^r2n+i is much smaller than the preceding term, while this is not the 
case for r < r^. This is most likely because the convergence is lost for r < r^; in the method of 
solution used. 

5. It is necessary to check as well that the results don't depend on the Tq. Let us take 
e = 0.58 and F = 7/5: when tq doubles, k^ decreases by 10% and decreases by 20%. The 
difference stems from the fact that the flow become slightly compressed, while moving from the 
sonic surface to tq , and, most importantly, its direction slightly changes. For smaller e, the ro 
dependence of the parameters of the region is even weaker. 

Note that the satisfactory accuracy of the results cannot be improved signiflcantly by aban- 
doning the assumptions made in favor of more accurate calculations. First, we used the flrst 
approximation for our calculations with a flnite e, and, second, we used its form at small radii 
where the flrsts approximation grows compared to the zeroth solution. We cannot specify any 
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Figure 2. Temperature versus coordinates in the plane of variables (r/rg, ^^o) 
for the parameters F = 7/5, £ = 0.03, Cqo = 0.00001. The light area at the center 
has the highest temperature. 



initial conditions near the sonic surface, since the behavior of formula (5) is asymptotic for 
r ^ 0. 

To make progress in this field requires numerically solving an exact system of two nonlinear 
partial differential equations, one first-order equation and one second-order equation. In this 
case, wc must reveal features that are many orders of magnitude smaller than the outer radius 
of integration or even the sonic surface, which can be achieved only by using adaptive methods 
with a variable step. In addition, we must numerically calculate the parameters of the separatrix 
surface (the boundary of the causally connected region) at which its smooth passage is realized, 
solve the time-dependent problem of establishment of this regime, show that it is stable, and 
find the percentage of cases, in which the accretion without shocks is possible. Such calculations 
are far beyond the scope of this paper. 



6. PHYSICAL AND OBSERVATIONAL SIGNIFICANCE OF THE RESULT 



Noteworthy among the features of the axisymmetric flow on a compact object observed so 
far is the tail convergence of streamhnes behind a moving object during accretion in a region 
of cold dense dust, where star formation takes place. 

The effect discussed in this paper has not yet been observed. The parameter F = 7/5 has a 
very narrow validity range. For example, it is realized at the initial stage of accretion onto a 
compact object from dense clouds of molecular hydrogen before its dissociation (the equations 
with a constant F cannot be used at all when the dissociation occurs). 

Let us consider in more detail the accretion of a hot gas from H II regions with F = 5/3 
far from the massive object. The value of F is the same on the sonic surface, since the gas 
density and temperature at infinity and on the sonic surface differ by a small factor (Bondi 
1952). In this case, the angular momentum of the gas relative to the stellar center in a stable 
regime is very low. Although the electrons become relativistic and the adiabatic index decreases 
sharply F = 13/9 starting from a certain distance (Shapiro and Teukolsky 1985), no region of 
anomalous compression is formed in this case. 

Nevertheless, let us calculate the energy release in the region of streamhne convergence for 
nonphysical constant values of F < 5/3. The calculation will yield an upper limit for actual 
astrophysical objects. For our estimation, we assume that all of the observed radiation in the 
derived temperature range is only thermal bremsstrahlung (Shapiro 1973). Its intensity per 
unit volume is defined by the non-relativistic formula A = cn^T^I"^ with the coefficient c that 
depends weakly on the number density n and temperature T . The intensity of the radiation 
for spherically symmetric accretion is then 



The characteristic radial size of the region in our case is Ar Tx/^i and the tangential size is 
~ T'xOxkx, where Ox defines the maximum angle at radius ro, starting from which fiow lines 
contracts in l/kx- The energy release in our region is 



For our estimation, we set ^-^ = 0.5. The ratio of the released energies can then be written as 



It follows from (14) and from the tables that the ratio of the released energies increases as the 
radius, at which the streamlines are focused, decreases. Let us calculate it for the minimum 
possible e when the focusing still takes place. We take the lower rows from the tables of param- 
eters for the region of anomalous compression for F = 4/3 and F = 7/5. We find from (14) that 
Lx/Lq reaches 100; i.e., when the focusing is present, the Bondi-Hoyle accretion efficiency can 
increase up to a factor of 10^. For physically meaningful value F = 7/5, the calculated ratio 
(14) actually specifies only an upper limit for Lx/Lq, since the gas dissociates rapidly and is 
subsequently ionized as it heats up, which ensures a transition to the regime with F = 5/3 and 
is accompanied by a significant decrease in Lx- However, the proper allowance for the change 
in adiabatic index in the axisymmetric case is beyond the scope of this paper. Although the ve- 
locity dispersion of compact objects is much higher than the speed of sound in the surrounding 
cloud of gas and dust, the described regime with e <1 can take place in certain cases at a low 
relative velocity of the gas cloud and the compact object. Such a change in energy release and 
accretion spectrum must be taken into account when analyzing experimental data to estimate 




L, 



X 




(14) 



the parameters of the black hole itself moving in regions of gas compression and its surrounding 
medium. 
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